The Maximal Amount of Gravitational Waves in the Curvaton Scenario 
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The curvaton scenario for the generation of the cosmological curvature perturbation on large scales 
represents an alternative to the standard slow-roll scenario of inflation in which the observed density 
perturbations are due to fluctuations of the inflaton field itself. Its basic assumption is that the initial 
curvature perturbation due to the inflaton field is negligible. This is attained by lowering the energy 
scale of inflation, thereby highly suppressing the amount of gravitational waves produced during 
inflation. We compute the power-spectrum of the gravitational waves generated at second order in 
perturbation theory by the curvaton (isocurvature) perturbations between the end of inflation and 
the curvaton decay. An interesting property of this contribution to the tensor perturbations is that 
it is directly proportional to the amount of non-Gaussianity predicted within the curvaton scenario. 
We show that the spectrum of gravitational waves may be in the range of future gravitational wave 
detectors. 



PACS numbers: 98.80.cq 

Inflation [l], Q has become the dominant paradigm for 
understanding the initial conditions for structure for- 
mation and for Cosmic Microwave Background (CMB) 
anisotropy. In the inflationary picture, primordial den- 
sity and gravity-wave fluctuations are created from quan- 
tum fluctuations "redshifted" out of the horizon during 
an early period of superluminal expansion of the universe, 
where they are "frozen" [3, 4]. These perturbations at the 
surface of last scattering are observable as temperature 
anisotropy in the CMB. The last and most impressive 
confirmation of the inflationary paradigm has been re- 
cently provided by the data of the Wilkinson Microwave 
Anisotropy Probe (WMAP) mission which has marked 
the beginning of the precision era of CMB measurements 
in space 

Despite the simplicity of the inflationary paradigm, 
the mechanism by which cosmological adiabatic pertur- 
bations are generated is not yet established. In the stan- 
dard slow-roll scenario associated to single-field models 
of inflation, the observed density perturbations are due 
to fluctuations of the inflaton field itself when it slowly 
rolls down along its potential. When inflation ends, the 
inflaton oscillates about the minimum of its potential 
V{4>) and decays, thereby reheating the universe. As a 
result of the fluctuations each region of the universe goes 
through the same history but at slightly different times. 
The final temperature anisotropics are caused by the fact 
that inflation lasts different amounts of time in different 
regions of the universe leading to adiabatic perturbations. 

An alternative to the standard scenario is represented 
by the curvaton mechanism [|| 0] where the final curva- 
ture perturbations are produced from an initial isocurva- 
ture perturbation associated to the quantum fluctuations 
of a light scalar field (other than the inflaton), the curva- 
ton, whose energy density is negligible during inflation. 
The curvaton isocurvature perturbations are transformed 
into adiabatic ones when the curvaton decays into radi- 
ation much after the end of inflation. 
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Contrary to the standard picture, the curvaton mech- 
anism exploits the fact that the total curvature pertur- 
bation (on uniform density hypersurfaces) £ can change 
on arbitrarily large scales due to a non-adiabatic pres- 
sure perturbation which may be present in a multi-fluid 
system. While the entropy perturbations evolve inde- 
pendently of the curvature perturbation on large scales, 
the evolution of the large-scale curvature is sourced by 
entropy perturbations. 

During inflation, the curvaton energy density is negli- 
gible and isocurvature perturbations with a flat spectrum 
are produced in the curvaton field a, 
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where if* is the value of the Hubble rate during infla- 
tion. After the end of inflation, the curvaton field oscil- 
lates during some radiation-dominated era, causing the 
ratio between its energy density and the radiation energy 
density to grow, thereby converting the initial isocurva- 
ture into curvature perturbation. The energy density p a 
will then be proportional to the square of the oscillation 
amplitude, and will scale like the inverse of the locally- 
defined comoving volume corresponding to matter dom- 
ination. On the spatially flat slicing, corresponding to 
uniform local expansion, its perturbation has a constant 
value Spa/pa ~ {6o~k/77* ) , where ef* is the value of the 
classical curvaton field during inflation. 

The curvature perturbation £ is supposed to be negligi- 
ble when the curvaton starts to oscillate, growing during 
some radiation-dominated era when p a j 1 p oc a, where a 
is the scale factor. After the curvaton decays £ becomes 
constant. In the approximation that the curvaton decays 
instantly it is then given by 
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where r = (p a /p)D and the subscript D denotes the 
epoch of decay. The corresponding spectrum is @ 

P <="GS:)- <3) 

It is nearly scale-invariant under the approximation that 
the curvaton field is effectively massless during inflation. 

The generation of gravity-wave fluctuations is another 
generic prediction of an accelerated de Sitter expansion 
of the universe. Gravitational waves, whose possible ob- 
servation might come from the detection of the -B-mode 
of polarization in the CMB anisotropy Q , may be viewed 
as ripples of spacetime around the background metric 

g M u = a 2 (t) (dr 2 - (<% + h l3 ) dx l dx j ) , (4) 

where r is the conformal time. The tensor hij is traceless 
and transverse and has two degrees polarizations, A = 
±. Since gravity- wave fluctuations are (nearly) frozen on 
superhorizon scales, a way of characterizing them is to 
compute their spectrum on scales larger than the horizon. 
During a de-Sitter stage characterized by the Hubble rate 
H* , the power-spectrum of gravity- wave modes generated 
during inflation is 

where M p = (SttGn)- 1 / 2 ~ 2.4 x 10 18 GeV is the Planck 
scale. Detection of the _B-mode of polarization in the 
CMB anisotropy requires H* > 10 12 GeV @. 

What about the expected amplitude of gravity-wave 
fluctuations in the curvaton scenario? 

The curvaton scenario liberates the inflaton from the 
responsibility of generating the cosmological curvature 
perturbation. Its basic assumption is therefore that the 
initial curvature perturbation due to the inflaton field 
is fully negligible. In the standard slow-roll inflationary 
models where the fluctuations of the inflaton field cj> are 
responsible for the curvature perturbations, the power- 
spectrum of the curvature perturbation is given by 

where ~ 1 is the spectral index and e = (<f> 2 /2M 2 H 2 ) 
is the standard slow-roll parameter. Requiring that the 
contribution (JB|) is much smaller than the value required 
to match the CMB anisotropy imposes H t -C 10 -5 M p . 
This implies that the curvaton scenario predicts the am- 
plitude of gravitational waves generated during inflation 
§5§ far too small to be detectable by future satellite ex- 
periments aimed at obser ving the B-mode of the CMB 
polarization (see however [13]). 



This is not the full story though. Gravitational waves 
are inevitably generated at second order i n p erturba- 
tion theory by the curvature perturbations [TJ LH, [H| . 
This scalar-induced contribution can be computed di- 
rectly from the observed density perturbations and gen- 
eral relativity and is, in this sense, independent of the cos- 
mological model for generating the perturbations. The 
generation of course takes place after the curvature per- 
turbation is generated. 

In the standard scenario, where the curvature pertur- 
bation is produced during inflation, the production of 
tensor modes occurs after inflation when the curvature 
perturbations re-enter the horizon. 

In the curvaton scenario, the production of tensor 
modes through the curvature perturbations may occur 
only after the curvaton decays, i.e. after the isocurva- 
ture perturbations get converted into curvature fluctu- 
ations. The energy density of gravitational waves (per 
logarithmic interval) is given by 

^Gw(k,T) = -^-^V h (k,T) (7) 

and the one generated by the curvature perturbations 
results to be of order flaw — 10 -20 , for those modes that 
re-entered the horizon when the universe was radiation 
dominated [II [LI. 

What we will be concerned about in this paper is the 
generation of tensor modes by the curvaton perturbations 
between the end of inflation and the time of curvaton 
decay. In other words, we are interested in the tensor 
modes generated at second order when the perturbations 
are still of the isocurvature nature. This contribution 
may be larger than the one created by the second-order 
curvature perturbations after the curvaton decay. In this 
sense, the spectrum of tensor modes computed in this 
paper corresponds to the maximal possible amount of 
gravity waves within the curvaton scenario. 

An interesting aspect is that the contribution to the 
tensor perturbations turns out to be directly proportional 
to the possibly large amount of Non-Gaussianity (NG) in 
the CMB anisotropics which is predicted within the cur- 
vaton scenario [3, LHj] . NG is usually parametrized in 
terms of the the nonlinear parameter /nl and the lat- 
ter is predicted to be of the order of 1/r in the curva- 
ton scenario; present-day data limit | /nl, I to be smaller 
than about 10 2 , that is r > 1CT 2 Therefore, the 

present observational bound on the level of NG in the 
CMB can already put an upper bound on the amount 
of tensor modes induced by the curvaton perturbations. 
This relic gravitational radiation may be particularly rel- 
evant in view of the realization that space-based laser in- 
terferometers, such as the Big Bang Observer (BBO) and 
the Deci-hertz Interferometer Gravitational wave Obser- 
vatory (DECIGO), operating in the frequency range be- 
tween ^0.1 Hz and 1 Hz may achieve the necessary sen- 
sitivity pjj. 

The equation for the second-order gravitational waves 
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before curvaton decay can be written (neglecting the 
first-order vector perturbations) as 



where = a' /a is the Hubble rate, the prime stands for 
differentiation with respect to the conformal time and 
k 2 = 8ttGn- If we define the Fourier transform of the 
tensor perturbations as follows 
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where the polarization tensors 
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where cos 9 = k • p. The second-order tensor modes are 
generated when the various modes k enter the horizon. 
Meanwhile, the production ends when the curvaton de- 
cays. We will assume in the following that the whole 
generation of tensor modes takes place in the radiation- 
dominated epoch. This assumption is motivated by re- 
quiring that the NG induced by the curvaton is sizeable, 
which requires the curvaton energy density not to dom- 
inate by the time of decay. Since the Hubble rate is 
given by Ti. = 1/r, a given mode k enters the horizon at 
Tk = 1/fc. Indicating by kjj the mode which enters the 
horizon at the time of the curvaton decay, kp = a(r£>)r, 
where T is the decay rate of the curvaton, we may write 
the evolution of the scale factor as 



4(k) = -^(kJe^kJ+Si^Qc)), 
er.(k) = ^Mk^-e^k^k)) 



(10) 



are expressed in terms of orthonormal basis vectors e and 
e orthogonal to k, the projector tensor in Eq. ([H]) reads 
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■e lk - x 4(k)e A * m (k). (11) 



In Fourier space, the equation of motion for the gravi- 
tational wave amplitude (for each polarization) then be- 
comes 

h£ + 2W4 + k 2 h k = S(k, r), 
5(k,r) = An 2 J ^0^e +lm (k)p lPm Sa p (r)6a^ p (±^ 

The solution to this equation can be easily found to be 

Mr) = 7~T f T dr'g k (r', r)a(r')5(k, r'), (13) 

where gk(r' , t) is the appropriate Green function either 
for a radiation- or a matter-dominated period. 

We split the perturbations of the curvaton field into a 
transfer function piece T CT (fc,r) and the primordial fluc- 
tuation (5(7k, 

6a k (r)=T (T (k,T)5a k (14) 
with the primordial power-spectrum defined by 
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The power-spectrum of the second-order gravitational 
waves becomes 



V h (k,T) 
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Trading the curvaton decay rate with the temperature at 
decay To, we obtain 



k 



D 



10" 



GeV 



Hz. 
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After the end of inflation, the zero mode a of the curvaton 
field starts oscillating at r m = (l/fcc^r/m) 1 / 2 , where m 
is the curvaton mass. Let us first consider those per- 
turbations which enter the horizon when the zero mode 
cf of the curvaton decay is already oscillating, that is 
k <(m/T)^ 2 k D . In this range of wavenumbers, one can 
show that the curvaton perturbations scale as the zero 
mode, Sa^r) ~ a ^ a~ 3 / 2 . This allows to write 
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where 7jd is the value of the curvaton zero mode at 
the time of decay. If we introduce the variables x = 
\k — p| jk and y — p/k and uses the radiation-dominated 
Green function gk(r', r) = sin(fc(r — r'))/fc, it is easy to 
realize that the main tensor mode production happens at 
horizon entry, that is at t ~ r^. Therefore, the power- 
spectrum (TIB]) computed at horizon entry is 



V h (k,T k ) ~ 10 2 /nl (-) I dy [ +J dx 



x y 2 ( (l+y 2 ~x 2 ) 2 ^ 2 
x 2 \ Ay 2 



P ( {kx)P c (ky), 



(20) 

where we have made use of the relations r = 
(k 2 /3T 2 )(m 2 a 2 J ) and /nl ~ In the curvaton sce- 

nario the resulting curvature perturbation is nearly scale- 
invariant and we can take Vq ~ (5 x 10~ 5 ) 2 . The remain- 
ing integrals in (|20p are dominated by the momenta for 
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which x ~ y ~ (m/T) 1 ^ 2 (kD/k). We finally obtain an 
energy density today of gravitational waves given by 

valid for k D <k< {m/Y) 1 / 2 k D . 

For those perturbations which enter the horizon before 
the zero mode of the curvaton field starts oscillating, that 
is k > (m/r) 1 / 2 fc£ ) , the scaling is is Sa^ ~ a -1 for < 
r < fcr^j, that is till fc > ma. For fcr^ < r < tq the 

scaling is is Jok ~ a -3 / 2 . Meanwhile, the zero mode 
a remains frozen till the mass of the curvaton becomes 
larger than the Hubble rate at r ~ r m . Repeating the 
steps leading to Eq. (j2"0|) , we find at horizon entry 

/ r \ f°° f 1+y 
V h (k,r k ) ~ 10 2 / 2 L (~J ^ dyj^ ^ dx 

- ^( i - (1+ 4,r a;2)2 ) 2 ^(Mn(^ 

(22) 

which corresponds to an energy density today of 

■taw* 10- (£)"(!), (23, 

valid for fc > (rn/T^/^kz). Finally, let us note that 
the redshifted gravitational wave at later times is al- 
ways larger than S/k 2 . Therefore the power-spectrum 
of gravity waves produced by the curvaton fluctuations 
is always bigger than the one generated by the second- 
order curvature perturbations which is inevitably gener- 
ated when the cosmological perturbations acquire their 
adiabatic nature. 

There is a simple physical motivation for the fact that 
the amount of gravity waves generated by the curvaton 
decay is enhanced by powers of 1/r with respect to the 



one produced by ordinary second-order curvature pertur- 
bations. Indeed, being the final adiabatic perturbations 
generated by the curvaton isocurvature perturbations, 
the smaller the amount of curvaton energy density at 
decay is, the larger the curvaton fluctuations have to be: 
<5<7k oc Ck/ r i see Eq. The isocurvature perturbations 
giving rise to gravity waves are therefore parametrically 
larger at horizon entry than the second-order curvature 
perturbations. 

From our findings we deduce that the amount of grav- 
itational waves in the perturbative regime T < m can be 
as large as Hqw — 10~ 15 maximazing the observation- 
ally allowed NG in the CMB anisotropics. This is quite 
intriguing since such a spectrum is at the range of, e.g., 
BBO and DECIGO interferometers. To be in the right 
frequency range, between 10 _1 and 1 Hz, one has to im- 
pose that the curvaton decays at To < 10 8 GeV. Further- 
more, the correlated BBO interferometer proposal claims 
a sensitivity down to J]qw — 10~ 17 . This would require 
r > 10~ 2 m in units of (10 2 //nl) ■ As a final remark, 
we point out that the expressions (f2Tj) and (|23[) are also 
applicable to the so-called modulated reheating scenario 
in which the curvature perturbations are due to the fluc- 
tuations of some light field parametrizing the inflaton 
decay rate [l8| . In such a case though the NG nonlinear 
parameter /nl does not exceed unity [TH - 

In summary, we have computed the maximal amount 
of tensor perturbations which may be generated within 
the curvaton scenario. It is directly proportional to 
the square of the nonlincarity parameter /nl which 
parametrizes the NG in CMB anisotropics. The present 
observational bound on the level of NG in the CMB data 
sets already an upper bound on the amount of tensor 
modes induced by the curvaton perturbations. On the 
other hand, a possible future detection of a large NG 
signal in CMB anisotropies would suggest a quantity of 
gravity waves at range of future gravitational wave de- 
tectors. 
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